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ABSTRACT 

We report multi-frequency Very Large Array observations of three massive 
star formation regions (MSFRs) containing radio continuum components that 
were identified as broad radio recombination hne (RRL) sources and hyper- 
compact (HC) Hll region candidates in our previous H92a and H76a; study: 
GlO.96+0.01 (component W), G28.20-0.04 (N), and G34.26-I-0.15 (B). An ad- 
ditional HC Hll region candidate, G45. 07-1-0. 13, known to have broad H66a and 
H76a lines, small size, high electron density and emission measure, was also in- 
cluded. We observed with high spatial resolution (0''9 to 2''3) the H53a, H66a, 
H76a, and H92a RRLs and the radio continuum at the corresponding wave- 
lengths (0.7 to 3.6 cm). The motivation for these observations was to obtain 
RRLs over a range of principal quantum states to look for signatures of pressure 
broadening and macroscopic velocity structure. We find that pressure broadening 
contributes significantly to the linewidths, but it is not the sole cause of the broad 
lines. We compare radio continuum and dust emission distributions and find a 
good correspondence. We also discuss maser emission and multi-wavelength ob- 
servations reported in the literature for these MSFRs. 
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1. Introduction 

Hypercompact (HC) Hll regions represent an earlier evolutionary state than compact 
and ultracompact (UC) Hll regions. These very small (< 0.03 pc), high electron density 
(ne ~ 10^~^ cm~^) and high emission measure (EM > 10^ pc cm^^) nebulae are ionized 
from within by O or B stars and usually are associated with maser emission. Many, but 
not all, have unusually broad radio recombination lines (RRLs) with full widths at half 
maximum intensity (FWHM) > 40 km s"^ (De Free at al. 1996, 1997; Sewilo et al. 2004a). 
HC H II regions have rising spectral indices a (where oc i/") from centimeter to millimeter 
wavelengths with typical values of a ~1, which suggests a range of densities and optical 
depths in the ionized gas. Broad RRLs (BRRLs) and intermediate sloped (—0.1 < a < +2) 
power-law spectral indices may be associated with the age of HC Hll regions, appearing 
only during a small fraction of the lifetime of the hypercompact phase. There are very dense 
and compact H ll regions, however, that do not have intermediate power-law spectral energy 
distributions or broad radio recombination hnes (e.g., G34.26-I-0.15 A, Avalos et al. 2006; 
Sewilo et al. 2004a). 

Although possibly related to the age, the physical reason that some HC Hll regions 
do not present BRRLs is poorly understood. The broad linewidths suggest highly dynamic 
internal structures (e.g., bipolar outflows, expansion, stellar winds, accretion, disk rotation, 
shocks) as well as electron impact (pressure) broadening. High spatial resolution observations 
over a wide range of frequencies are required to distinguish the relative contributions from 
pressure broadening and bulk motion of the gas. Transitions with high principal quantum 
numbers, n, are significantly affected by pressure broadening, while low n transitions are 
mostly free of pressure broadening and are sensitive to large-scale motions. With high spatial 
resolution at a range of frequencies, it is possible to distinguish bulk motions from pressure 
broadening. Observations have shown that both pressure broadening and gas dynamics 
contribute to hne broadening (Sewilo et al. 2008; Keto et al. 2008). 

Theoretical progress has been made towards understanding the nature of HC H II re- 
gions. The model of hierarchical clumping of the nebular gas proposed by Ignace & Church- 
well (2004) was able to reproduce the intermediate values of spectral indices seen in HC 
H II regions. Tan & McKee (2003) and Keto (2003) proposed models that predict both con- 
finement of the ionized gas in HC H lis and the BRRL emission. The photoevaporating disk 
model, originally proposed by HoUenbach et al. (1994) and Lizano et al. (1996) to explain 
the most compact UC Hll regions, has been recognized to be applicable to HC Hll regions 
(Lugo et al. 2004). 
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Here we report high spatial resolution, multi-frequency Very Large Array^ (VLA) ob- 
servations of four candidate HC Hll regions with the goal of determining their physical prop- 
erties. Three sources, GlO.96+0.01, G28.20-0.04, and G34.26+0.15, were selected based on 
an earlier H92q; and H76q; study (Sewilo et al. 2004a). We also observed the HC Hll can- 
didate G45.07+0.13 NE, which is known to have broad H66a and H76q; lines (Garay et al. 
1985, 1986), small size, and high electron density (i.e., high emission measure). Additional 
RRLs were observed to search for systematic broadening with increasing principal quantum 
number. A strong correlation of line width with principal quantum number would indicate 
pressure broadening as the main mechanism for line broadening, as opposed to large-scale 
motions (rotation, outflows, infall, etc.). The rms electron temperature and the rms electron 
density of the sources can be determined for each transition using the assumption of local 
thermodynamic equilibrium (LTE). RRLs also provide information on the velocity structure 
of the sources. 

We report detections of the H53q;, H66q;, H76q;, and H92q; RRLs and the 0.7, 1.3, 2.0, 
and 3.6 cm radio continuum emission with resolutions of 1''5, 0''9, 1''2, and 2''3, respectively, 
toward the MSFRs GlO.96+0.01, G28.20-0.04, G34.26+0.15, and G45.07+0.13 (hereafter 
G10.96, G28.20, G34.26, and G45.07). The observations are described in Section 2. In 
Section 3, we derive physical properties of the nebulae based on continuum observations, 
and investigate the change of line width with principal quantum number. In Section 4 we 
investigate the dust emission toward these sources using mid- infrared Spitzer Space Telescope 
GLIMPSE data. We compare the mid-infrared GLIMPSE and radio continuum images with 
previous maser emission and multiwavelength observations. We discuss the general results 
in Section 5 and summarize the main conclusions in Secion 6. 



2. The Data 

2.1. VLA Observations: Continuum and Radio Recombination Lines 

The VLA was used to observe H53q;, H66q;, H76q;, and H92q; RRLs toward G10.96, 
G28.20, G34.26, and G45.07. Table 1 gives the instrumental parameters of the observations. 
Table 2 summarizes the transitions that we observed toward each source (Sewilo et al. 2004a, 
and this paper) and the adopted distances. 

The first channel of VLA spectral line data ("channel 0") contains the average of the 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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central 75% of the observing bandwidth (generated by vector averaging channels; ?). The 
broad-band uv data from "channel 0" were edited and calibrated using the Astronomical 
Image Processing System (AIPS) software package, and the resulting calibration solutions 
were applied to the multi-channel spectral-line data. All sources were self-calibrated to in- 
crease the dynamic range of the images. The continuum level was determined from a linear 
fit to the line-free channels. Continuum subtraction was done in the uv plane with the AIPS 
task UVLSF. A line-free continuum data set and a 31 channel cube of recombination line 
emission were constructed. Both data sets were Fourier inverted using weighting intermedi- 
ate between natural and uniform (Robust=0 in AIPS task IMAGR) and CLEANed in the 
standard manner. The images were corrected for primary beam attenuation using the AIPS 
task PBCOR. Integration times and continuum noise values of the final images are listed in 
Table 3. 

Source positions, deconvolved angular diameters, position angles, and total flux densities 
(Table 3), were determined by elliptical Gaussian fitting using the least-squares AIPS routine 
JMFIT. Continuum-subtracted line images were used to generate the spatially integrated 
H53q;, H66q;, H76q;, and H92q; line profiles. For each line, the integration was performed 
approximately over the area enclosed by the 50% contour with respect to the peak value. The 
subsequent spectral analysis of the RRLs was done using the IDL Programming Language 
^. The spatially integrated hne profiles were fitted by a single Gaussian. The results for 
individual sources are presented graphically in Figs. 1 - 6 for the H53q;, H66q;, H76q;, and 
H92q; RRLs. The continuum images are shown in the same figures; the projected linear 
scales are based on the kinematic distances listed in Table 2. 



2.2. Archival Data: VLA and Spitzer Space Telescope 

In Sewilo et al. (2008), we reported the results of high- resolution (0''15) observations of 
G28.20 N that revealed the shell-like morphology of this HC Hll region. To investigate the 
morphology of G28.20 N at even higher resolution, here we use VLA K-band A-array archival 
data (Program ID AZ168) to create the high resolution (0''09) continuum image shown in 
Fig. 4. These observations were taken on 2006 March 14, with a 3.125 MHz bandwidth of 63 
channels of 48.828 kHz each. We used line-free channels of these spectral line data to make 
a self-calibrated continuum image with an rms noise of 85 fiJj beam~^. The measured 1.3 
cm flux density of G28.20 N is 630±63 mJy The synthesized beam is 0''10x0''08 at position 
angle of — 0?5. 



^www.ittvis.com/idl 
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We also present H2O, OH, and CH3OH maser positions relative to the radio continuum 
for each HC H ll region based on high resolution (interferometric) data from the literature. 
The relative positional accuracy of 22 GHz H2O masers detected toward G45.07 NE is O'.'l 
(Hofner & Churchwell 1996). The 1665 MHz and 1667 MHz OH masers observed in G10.96 
W and G28.20 N have an absolute positional accuracy of ~0''3 in right ascension and ~0''5 
in declination (Argon et al. 2000). The absolute positions of the 6.67 GHz CH3OH masers 
in G10.96 W and G28.20 N are accurate to ~1" (Walsh et al. 1998), while the positional 
uncertainty of the 44 GHz CH3OH masers in G45.07 NE is 0% (Kurtz et al. 2004). The 
relative positional accuracy of 22 GHz H2O masers detected toward G45.07 NE is 0''1 (Hofner 
& Churchwell 1996). The distribution of masers is shown in Figs, lb, 4, and 6 for G10.96 
W, G28.20 N, and G45.07 NE, respectively. 

We used VLA K-band spectral-hne archival data for G28.20 N HC Hll region (Program 
ID AP465) to supplement the single-dish water maser data available for this source in the 
literature. The observations were taken on 2008 April 4 with the VLA in C configuration. A 
3.125 MHz bandwidth with 255 channels of 12.2 kHz each was used. We followed standard 
VLA data reduction procedures. The synthesized beam is l''19x0''94 at position angle of 
— 12?8. The typical rms of a line-free channel is 22 mJy beam"^. We determined the positions 
and peak flux densities of the masers from Gaussian flts to the peak channel of each maser. 
The H2O maser positions and observed parameters are listed in Table 4. The uncertainty of 
the H2O maser positions in G28.20 N is O^'l. 

We also investigate dust emission for our sample using data from the Spitzer Space 
Telescope GLIMPSE survey ("Galactic Legacy Infrared Mid-Plane Survey Extraordinaire", 
Benjamin et al. 2003; Churchwell et al. 2009). The inner Galactic plane was imaged at 
3.6, 4.5, 5.8, and 8.0 //m using the IRAC camera (Fazio et al. 2004) on the Spitzer Space 
Telescope. The IRAC spatial resolution ranges from 1''5 to 1''9 between the 3.6 and 8.0 //m 
bands. The IRAC image mosaics were created by the GLIMPSE pipehne. The absolute 
accuracy of the GLIMPSE point source positions is 0''3 (Meade et al. 2009). 

3. Results 

3.1. Physical Properties of the Sources 

Physical properties of the sources were derived from the observed continuum param- 
eters reported in Table 3 using the equations of Mezgcr & Henderson (1967) and Panagia 
& Walmsley (1978) for spherical, optically thin, homogeneous, and ionization-bounded Hll 
regions. Panagia & Walmsley (1978) presented modifled versions of the Mezger & Henderson 
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(1967) equations for geometry-dependent physical parameters such as the rms electron den- 
sity {ne^rms), the emission measure (EM), and the mass of ionized gas (Mhh). The excitation 
parameter (U) and the Lyman continuum photon flux required to maintain ionization of the 
nebula (N^) do not depend on the source size or its morphology and were derived using the 
formulas from Mezger & Henderson (1967). The continuum optical depth (tc) was calculated 
from the Oster (1961) formula in a form presented by Mezger & Henderson (1967), using EM 
derived based on the results of Panagia & Walmsley (1978). The derived physical parameters 
for each source are presented in Table 5. The hnear diameter of a sphere hsted in Table 5 
is the diameter of the spherical model for the source, rather than the apparent Gaussian 
half-power width of the source reported in Table 3. An optically thin homogeneous sphere 
has the same peak EM and total flux density as the observed source. The calculations take 
into account the dependence on the beam size as described in Panagia & Walmsley (1978). 

The values were used to estimate a lower hmit to the spectral type of the excit- 
ing star given in the last column of Table 5. The classification is based on the non-LTE, 
spherically expanding (owing to the stellar wind), line-blanketed (to account for the effects 
of metals on the atmospheric structure and emerging ionizing flux) model atmospheres de- 
rived by Smith et al. (2002) under the assumption that a single star ionizes a dust-free, 
ionization bounded H ll region. The assignment of a given spectral type with depends 
strongly on assignment of effective temperature with spectral type and how accurately stel- 
lar atmosphere models reflect reahty. For these reasons, spectral types hsted in Table 5 are 
approximations. The stellar spectral type corresponds to the value of radio flux density at 
the highest frequency observed. 

A second set of physical properties of the nebulae was derived based on their peak 
brightness temperature; these are "peak physical properties" (Wood & Churchwell 1989). 
These quantities depend only on the peak brightness, not the measured diameter or mor- 
phology. Thus, they can be derived for all the sources from our sample including the blended 
components B and C in the G34.26 complex and the over-resolved component E of G10.96. 
The brightness temperature Tb is given by the Rayleigh- Jeans approximation. 

The peak optical depth was estimated from Tb — Te{l — e~'^'') assuming a kinetic 
temperature of Te = 10^ K that is constant over the synthesized beam. The peak emission 
measure was calculated using the optical depth equation from Mezger & Henderson (1967). 
The relation =< >V2= (EM/AsY^^ was used to calculate the rms electron density 
(^e,rms) avcragcd over the beam. The quantity As is the source size along the line of sight 
through the peak position. We assume the sources in our sample have spherical geometry 
based on their projected circular symmetry; for As we adopt their geometrical mean angular 
diameters (half-power widths, or "source size" in Table 3) assuming that the extent of the 
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source along the line of sight is the same as its extent in the plane of the sky. The peak 
physical properties are given in Table 6. 



3.2. Radio Recombination Line Analysis 

The Gaussian fit parameters of the hne profiles are given in Table 7. Vlsr is the 
central line velocity relative to the local standard of rest, FWHM is the hne full width at 
half maximum, Sl is the flux density in excess of the continuum at line center, Sc is the 
continuum flux density, |^ and J ^ dv are the line-to-continuum ratios at line center and 
integrated over the line profile, respectively. The continuum fiux densities are determined 
by integrating over the same areas of the source used to generate the line profiles; thus 
they may be less than the fiux densities given in Table 3 (integrated over the entire source). 
The uncertainties given in Table 7 are the formal la (68.3% confidence level) errors of the 
nonlinear least-squares fit. These uncertainties are lower hmits because they do not account 
for systematic errors. 

The measured line widths given in Table 7 are broadened because of the finite velocity 
resolution of the observations. To determine the intrinsic line widths, the instrumental 
spectral resolution was deconvolved from the observed lines. The convolution function in 
the spectral frequency domain was assumed to be a Gaussian function with FWHM of 
1.2 times the channel width (or twice the channel width if the spectrum was Hanning- 
smoothed off-line) and deconvolved from the observed line widths. The deconvolved line 
widths (FWRMdeconv) are given in Table 8. 

We calculated the LTE electron temperatures T* for all the sources with reliable RRL 
measurements using the formula from Rohlfs & Wilson (1986): 

6.985 x 10=^ {SL/Sc)p,\k 1 



^" a(l/,Te) ' FWHMoecaav 1 + Y+. 

where v is in GHz, FWHMdecanv is in km s-^ and T* in K. F+ = N{He+)/N{H+) is the 
ratio of ionized He to H atoms and is assumed to be 0.1. The factor a(i/, T^) is approximately 
1 (Altenhoff et al. 1960). This equation was derived under several assumptions about Hll 
regions: (1) the structure is plane-parallel, homogeneous and isothermal; (2) all optical 
depths are small: [t; -|- rd <^ 1 and Tc <^ 1; and (3) the lines can be treated in the LTE 
approximation. 

However, several of the observed sources have relatively high optical depths {tc from 
0.02 to 0.2), causing attenuation of the line-to-continuum ratio. In these cases, the LTE 
electron temperatures have been calculated using the line-to-continuum ratios corrected for 
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opacity rather than the measured values: 

{SL/Sc)corr = {SL/Sc)peak ' (e^= - 1)/Tc. (2) 

The corrected hne-to-continuum ratios and the LTE electron temperatures are listed in Table 
8. The uncertainties in T* were estimated by a linear propagation of errors. 

3.3. Spitzer /IR AC Mid-Infrared Emission 

High-spatial resolution mosaics from GLIMPSE (Churchwell et al. 2009) are a unique 
tool for the study of massive star formation regions. The IRAC 3.6, 5.8, and 8.0 iim bands 
contain strong polycyclic aromatic hydrocarbon (PAH) features that probe conditions in 
photodissociation regions (PDRs) that separate the ionized gas in H ll regions from the sur- 
rounding molecular gas (e.g., Tielcns ct al. 1993). PAHs are large organic molecules contain- 
ing >50 C atoms and are very abundant in the ISM (e.g., Sellgren 1984; Allamandola et al. 
1989). PAHs are excited by sub-Lyman ultraviolet photons but are destroyed in Hll regions 
by extreme ultraviolet (EUV) radiation (Volt 1992; Povich et al. 2007). Nevertheless, Hll 
regions can be traced by the IRAC 4.5 /xm band, which does not contain PAH features. 
The 4.5 yum emission of Hll regions is a combination of thermal emission from small grains, 
scattered star light, and the Bra recombination line. The 3.6 and 4.5 /im bands are sensi- 
tive stellar tracers. The stellar contribution at 5.8 and 8.0 //m is relatively low because the 
photospheric emission is on the Rayleigh- Jeans tail of the energy distribution and declines 
rapidly with increasing wavelength. 

We present IRAC images in Figs. 7-10. All the images are presented in Galactic coor- 
dinates and the intensity is scaled logarithmically. The positions of radio sources are marked 

or the radio contours are overlaid on the images for reference. Equatorial coordinates are 
indicated to facilitate comparison with the radio images shown in Figs. 2-6. The GLIMPSE 
Point Source Catalog (green circles) and Archive^ sources (red circles) are overlaid on the 
3.6 fim images which show the stellar distribution in each field. UC and HC Hll regions are 
usually bright in the IRAC bands; thus they are often saturated in the IRAC images. G28.20 
N is saturated in the 8.0 and 4.5 iim bands, G34.26 is saturated in 8.0 and 3.6 /xm bands, 
G45.07 is saturated in all bands. G10.96 and G28.20 S are not saturated in any band. The 
fact that the sources are saturated does not affect the conclusions we draw in this paper. 



^The GLIMPSE Catalog is more reliable but less complete than the Archive. The GLIMPSE Data 
Delivery documents are available at http://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/ 
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The bright IRAC sources suffer from the "electronic bandwidth effect"^. This instru- 
mental effect is exhibited only in the 8.0 and 5.8 /im IRAC bands for saturated or bright 
sources; it produces artifacts appearing as a sequence of sources of decreasing intensity. This 
effect is seen in the 8.0 fim image of G28.20 N, and 5.8 /im and 8.0 fim. images of G34.26 
and G45.07 NE. As an example, the source appearing to be blended with G28.20 N slightly 
to the south in the 8.0 /xm image of Figure 8 is not the S component, but rather is an image 
artifact produced by this effect. 

A discussion of the IR emission from individual sources is presented in Section 4. 

4. Comments on Individual Sources 

In this section, we present the results of our RRL and mid-IR study for the individual 
sources. We also summarize the previous work on radio, IR, maser, and molecular line 
(where available) emission from each source. 

4.1. GlO.96+0.01 

The G10.96 MSFR consists of two compact radio continuum components G10.96 E 
and G10.96 W surrounded by more diffuse ionized gas. They have been detected in radio 
continuum between 1.4 GHz and 8.6 GHz (Becker et al. 1994; Walsh et al. 1998; Paper I). In 
Paper 1, we reported the detection of a broad H92« line with a FWHM of 43.8 ±1.5 km s~^ 
toward the W component. The E component has a lincwidth of 33.6 ± 2.8 km s^^, typical 
for UC Hll regions. Based on the 3.6 cm and H92a observations, we identified G10.96 W as 
an HCll region candidate (Paper I). 

The deconvolved linewidths of the H76q; and H66q; lines from G10.96 W reported in 
Table 8 are 28.4 ± 1.5 and 28.0 ± 0.6 km s^-*^, respectively (see Figure 1). These results 
indicate that the H92q; line is strongly affected by pressure broadening (see the discussion in 
Section 5.1). For G10.96 E, both H76q; and H66q; hues are narrow with FWHMs of 19.4 ± 6.5 
and 15.0 ± 4.9 km s~^, respectively (see Figure 2). These hues are narrower than their H92q; 
counterpart by over 10 km s~^, even if the 30% uncertainty in the FWHMs of the higher 
frequency hues is accounted for. The difference in FWHMs over a range of RRL transitions 
in a relatively low density gas indicates that other line broadening mechanisms may be at 
work in addition to the thermal and turbulent broadening and a small contribution from 



^See http://ssc.spitzer.caltech.edu/irac/features.html#2C 
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pressure broadening (see Section 5.1). Most of the diffuse 3.6 cm radio continuum emission 
observed near component E (Paper I; HPBW ~ 2''3) is resolved out in our 1.3 and 2 cm 
images (see Figure 2). 

Several maser species have been detected toward G10.96. Methanol (CH3OH) 6.67 GHz 
masers were first reported by Schutte et al. (1993) and later confirmed by Walsh et al. 
(1997), both based on single-dish observations. The high resolution interferometric survey 
of Walsh et al. (1998) revealed two CH3OH maser spots in the W component, near the con- 
tinuum emission peak (see Fig. lb). Water masers and CS(2-1) emission arc reported by 
CodcIIa et al. (1995) and Bronfman et al. (1996), respectively, based on single-dish observa- 
tions. 

In Figure 7 we present 3.6 and 8.0 images of the G10.96 MSFR. These images indicate 
a close correspondence between the distribution and morphology of dust and ionized gas. 
The diffuse dust emission from G10.96 E is filamentary, with a size of ~60"x40" (~4.1x2.7 
pc at a distance of 14 kpc; Table 2) in RA and Dec, respectively. This filamentary structure 
is more prominent at 24 /xm (MIPSGAL survey; Garey et al. 2009), but is also noticeable 
at shorter wavelengths. In the 3.6 /im image a possible central ionizing source (high-mass 
(proto)star) is prominent in the W component. The 3.6, 4.5, 5.8, and 8.0 /im flux densities 
of the IRAC source (GLIMPSE Archive ID: SSTGLMA GOlO.9584+00.0219) are 9.7 ± 0.9 
mJy, 58.7 ± 7.5 mJy 174.2 ± 12.1 mJy and 395.3 ± 30.7 mJy respectively. 

The IRAG images also reveal a cluster of sources around the radio emission peak in 
component E (see Figure 7). Some of these objects may be ionization sources; however, 
wc can not rule out that some may be background objects. The multiplicity of ionizing 

(proto)stars may give rise to the multiple-peak morphology of the UC Hll region shown 
in Figure 2. The mid-IR source located near the 3.6 cm peak of G10.96 E in Figure 7 
corresponds to one of the two peaks (E2 to the south-west) of the 2 cm emission (Figure 2b). 
The GLIMPSE Archive provides the 3.6, 4.5, and 5.8 /xm flux densities for this source 
(SSTGLMA GOlO.9646+00.0098) of 8.2 ± 2.7 mJy, 13.1 ± 4.0 mJy, and 58.9 ± 7.6 mJy, 
respectively. 

4.2. G28.20-0.04 

The G28.20— 0.04 MSFR hosts two radio continuum components: N and S. Both com- 
ponents were detected in previous radio continuum (0.7 cm - 6 cm, e.g. Kurtz et al. 1994; 
Walsh et al. 1998; Sewilo et al. 2004a, 2008) and RRL (H92q; and H53q;; Sewilo et al. 2004a, 
2008) observations. The N component has been identifled as the HG H ll region candidate 
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by Sewilo et al. (2004a) based on the 3.6 cm continuum and H92ci; RRL observations. The 
H92q; Unes reported in Paper I have FWHMs of 74 ± 3 and 35 ± 1 km for components 
N and S, respectively. Higher-frequency RRLs from G28.20 N (H53q;, Sewilo et al. 2008; 
H30q;, Keto et al. 2008) are significantly narrower than the H92q! line, with line widths of 
39.8 ± 1.7 km s~^ and 20.9 ± 0.6 km s~^, respectively. 

In Fig. 3, we present a 2 cm and a new 7 mm continuum image of G28.20, along with 
the corresponding H76a and H53a line profiles. Both RRLs are narrower than the H92« 
line reported in Paper I for both the N and S components. The H76q; and H53« FWHMs 
arc 55 ± 3 and 38 ± 1 km for component N, and 24 ± 4 km and 19 ± 2 s^^ for 
component S, respectively. The significant increase in line width with increasing quantum 
number for G28.20 N indicates that pressure broadening is an important (possibly dominant) 
line broadening mechanism for the H92q; line detected from this source. 

Although not evident in Figure 3, we detect marginal evidence for 7 mm continuum 
emission at the south-east border of G28.20 N. The emission is slightly more evident in the 
6 cm CORNISH image (C. Purcell, personal communication). We do not report parameters 
for this marginal detection, but note the possible presence of a non-thermal source at this 
position. In addition, our continuum images — at both 7 mm and 2 cm — show evidence 
for a new source, with a rising spectral index, just to the cast of G28.20 S. We refer to this 
source as G28.20 SI (see Figure 3a), and report approximate values for its parameters in 
Table 3. 

However, Sewilo et al. (2008) showed that gas dynamics also contributes to hne broad- 
ening in G28.20 N. The morphology and velocity structure of G28.20 N was resolved in their 
high resolution (~0''15) 7 mm continuum and H53q; VLA observations. The source appears 
as a shell-like structure with an inner radius of 1100 AU and an outer radius of 2500 AU. 
They detected a velocity gradient of 10^ km s^^ pc^^ along the minor axis of the continuum 
source, from which they inferred a rotating torus around an ~30 M0 central object. Based 
on line broadening and splitting along the major axis, they also suggested the presence of 
an outflow. 



4.2.1. Molecular Gas Kinematics 

G28.20 N has also been the target of high resolution molecular line observations by 
Qin et al. (2008), who proposed the existence of a hot core in G28.20 N. They used the 
SMA to image transitions in CH3CN, CO, ^^CO, SO2, OCS, and CH3OH with 1" resolution. 
These authors derived a CH3CN rotation temperature of 308 ± 22 K, inferred an H2 density 
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>10^ cm"^, and an age of G28.20 N of 1.5 x 10^ years. These properties, together with the 
small size (<0.04 pc) fulfill the standard hot core criteria (size <0.1 pc, density ^10^ cm~^, 
and temperature >100 K; Kurtz et al. 2000). Spectra of SO2, OCS, and CH3OH have similar 
radial velocities and line widths as CH3CN; thus they likely originate in the hot molecular 
core. Qin et al. (2008) argue that ^^CO emission originates outside of the hot molecular 
core, since the excitation temperature and H2 density are lower than those derived from the 
CH3CN hne. The ^^CO spectrum shows blueshifted absorption and redshifted emission with 
the absorption and emission peaks ahgned with the source rotational axis (i.e., the SE-NW 
direction; Sewilo et al. 2008). Based on the ^^CO distribution and hne profile, Qin et al. 

(2008) concluded that ^^CO traces an outflow from the molecular core over a hnear scale of 
~0.1 pc. 

The OCS position-velocity diagram for the NE-SW direction across the continuum 
source G28.20 N shows two velocity peaks (95 and 97 km s~^) separated by 0('2. This velocity 
difference corresponds to a NE-SW velocity gradient of ~400 km s~^ pc^^ (Qin ct al. 2008). 
A similar velocity pattern was observed in the K=5 component of the CH3CN (J=12-ll) 
transition detected toward G28.20 N (S.-L. Qin, private communication). The infall and 
outflow were previously reported by SoUins et al. (2005), based on VLA observations of 
ammonia (HPBW from ~0^'3 to ~3"). 

Based on the same SMA data previously used by Qin ct al. (2008), Klaassen et al. 

(2009) reported a detection of the SO2 and OCS velocity gradient of ~150 km s^^ pc^^ in 
a roughly N-S direction (as opposed to the NE-SW direction suggested by Qin et al. 2008). 
Klaassen et al. (2009) also suggested the presence of a CO outflow in the NE-SW direction. 
Qin et al. (2008) found that the CO data are inconclusive due to the contamination from 
the infall; as mentioned above, they suggested an outflow in the SE-NW direction based 
on the ^^CO data. The detection of the H53q; velocity gradient in the NE-SW direction by 
Sewilo et al. (2008) supports the Qin et al. (2008) interpretation. 

4.2.2. Maser Emission 

Numerous masers have been detected toward the G28.20 MSFR: 1665 and 1667 MHz 
OH (Argon et al. 2000), 22 GHz H2O (Han et al. 1998; Kurtz & Hofner 2005; this paper), and 
6.67 GHz CH3OH (Walsh et al. 1998). G28.20 was also observed with arcminute resolution 
in the excited OH transitions at 6035 MHz (Caswell & Vaile 1995b; Caswell 2003) and 
6030 MHz (Caswell 2003). The OH, H2O, and CH3OH maser emission is associated with 
component N (see Fig. 4). 
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4-2.3. Dust Emission 

Spitzer IRAC images of G28.20 at 8.0 and 3.6 /iin are shown in Fig. 8. The 8.0 /xm 
image shows that G28.20 is located toward an elongated region of high extinction (infrared 
dark cloud, IRDC) stretching over at least 2' (~3.3 pc at 5.7 kpc) in a roughly NE-SW 
direction. This IRDC is probably the natal molecular cloud from which the protostars in 
G28.20 are forming. Figure 8b shows the 3.6 //m image of G28.20 MSFR. The contours of the 
2 cm radio emission indicate the position of the N and S components. The mid-IR sources 
from the GLIMPSE Catalog/Archive are indicated by circles. A bright source, detected in 
Ks-band and all four IRAC bands, coincides with G28.20 N. A much weaker source, detected 
at 3.6 and 4.5 /xm, coincides with G28.20 S. Midway between the two radio components, a 
third IR source is a near-IR (JHKs) source (also detected at 3.6 fim). This source does not 
have a radio counterpart and may be a foreground object or a less deeply embedded object. 

Both G28.20 N and G28.20 S have counterparts in the GLIMPSE Archive: SSTGLMA 
G028.2003-00.0493 and SSTGLMAG028. 1984-00.0501, respectively. The mid-IR source cor- 
responding to G28.20 N is significantly brighter than the sources that matches G28.20 S. The 
3.6 and 4.5 /xm flux densities are 358.4 ± 18.9 mJy and 2123.0 ± 99.3 mJy for component N, 
and 2.1 ± 0.3 mJy and 8.9 ± 1.2 mJy for component S. The 5.8 //m flux density of G28.20 
N is 3892.0 ± 98.6 mJy. 

4.3. G34.26+0.15 

Four radio continuum components (A-D) and two infrared sources (E and F, Camp- 
bell et al. 2000) have been identified in the G34.26 MSFR. Component C is a prototypical 
cometary UC H ll region with a compact head toward the east and a diffuse tail to the west 
(e.g., Reid & Ho 1985; Garay et al. 1986; Sewilo et al. 2004a). Components A and B, located 
to the east of component C, are HC Hll region candidates. A more evolved component (D), 
located south-east of components A-C, is thought to be an Hll shell with a diameter ~1', 
expanding within an ambient density gradient (Fey et al. 1994). 

We detected H53q; emission toward G34.26 components A, B, and C. The fine parame- 
ters are presented in Tables 7 and 8; line profiles and the 7 mm continuum image are shown 
in Figure 5. 

The H76q; (Paper I) and H53q; (this paper) fines from G34.26 A are narrow, with 
FWHMs of ~22 km s""^, consistent with a combination of thermal and turbulent broadening. 
The H53q; linewidth of component C is 50.2 ± 6.4 km s"^, in good agreement with the H76q; 
linewidth of 55.5 ± 1.9 km s~^ reported by Garay et al. (1986). These authors demonstrate 
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that this relatively broad linewidth probably results from microturbulent motions having 
a velocity of ~30 km s~^. Component B has a broad H76q; line of FWHM 48.4±4.4 km 
(Paper I); however the H53q; line is ~15 km s~^ narrower (see Table 8) indicating that 
the H76q; line has a substantial contribution from pressure broadening. Also notable is the 
higher H53q; line center velocity for component B; the velocity is more than 15 km s~^ higher 
than for either component A or C. It is unclear if this velocity difference is important to 
understand the nature of component B. Hatchell et al. (2001) considered component B as a 
strong candidate for the driving source of the SiO outflow; if this is the case, there may be 
a dynamic explanation for the higher H53q; velocity. 

If we assume no contribution to the component A and B line widths from large scale 
motions and a non-thermal (turbulent) contribution of 5 km s~^ (FWHM; Cesaroni et al. 
1991, Plume et al. 1997, Sewilo et al. 2008), then we can determine an upper hmit for pressure 
broadening from the H76q; line. The electron temperatures for the A and B components are 
5700 K and 5300 K. respectively (Paper I). Thus, the thermal plus turbulent contribution to 
the line width is ~17 and ~16 km s~^ for components A and B, respectively. The resulting 
upper limits for the pressure broadening contribution to the II76q; linewidth (calculated as 
described in Sewilo et al. 2008) are ~10 km s~^ and ~43 km s~^ for the A and B components, 
respectively. 

Components A and B meet the criteria for HC H II regions (Paper I; Avalos et al. 2006; 
this paper). The diameters, electron densities and emission measures derived in Paper I 
based on high-resolution (~ 0''5) 2 cm observations are 0.008 pc and 0.006 pc (assuming a 
distance of 3.7 kpc); 1.4 x 10^ cm~^ and 2.2 x 10^ cm"^; 2.2 x 10^ pc cm"*^ and 4.3 x 10^ 
pc cm~^, respectively. However, both components A and B have a relatively low turnover 
frequency of ~ 10 GHz (Avalos et al. 2006). The physical parameters derived in this paper 
(see Table 5) for component A agree with the Paper I results. No physical parameters for 
component B can be derived because this source is partially blended with component C at 
the resolution of 1'.'5. 

Avalos et al. (2009) resolved G34.26 A and B with 0''05 resolution VLA observations and 
detected possible limb brightening in both sources. They used spherical shell models with 

power-law density gradients {n^, oc r~") to explain their observations. They showed that the 
7 mm intensity profiles, radio continuum spectra, and angular sizes of both G34.26 A and B 
can be reproduced by a shell of inner and outer radii ~400 AU and ~1000 AU, respectively, 
and a density gradient a ~ 0.3-1.0. This physical picture is quite similar to what has been 
proposed for G24.78+0.08 Al (Beltran et al. 2007) and G28.20-0.04 N (Sewilo et al. 2008). 
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4.3.1. Maser Emission 

Mascr emission is detected in G34.26 including H2O (Hofner & Churchwell 1996), OH 
at 1665 and 1667 MHz (Garay et al. 1985; Gaume & Mutel 1987, Gasiprong et al. 2002), 
OH at 6 GHz (Caswell k Vaile 1995b), and CH3OH at 6.67 GHz (Menten 1991) and 44 
GHz (Kurtz et al. 2004). The H2O maser spots are distributed ahead of and projected onto 
the cometary arc of component C, while most of the OH masers are more tightly confined 
to a parabohc arc along the eastern edge of component C (e.g., Gasiprong et al. 2002). A 
cluster of OH 1665 MHz masers is associated with component B. No OH, H2O, or CH3OH 
maser emission is detected toward component A. Several of the authors cited above provide 
detailed discussions of the relation between the masers and other physical phenomena within 
the MSFR. 

High angular resolution molecular line observations (Keto et al. 1992; Garay & Rodriguez 
1990; Heaton et al. 1989; Hofner & Churchwell 1996; Mookerjea et al. 2007) indicate that 
the molecular gas lies roughly between components A, B, and C (hot core). In particular, 
there is no evidence for localized molecular peaks coincident with components A or B. 

4.3.2. Dust Emission 

Mid-IR emission from components A, C, and D of G34.26 was detected at wavelengths 
from 8 to 20.6 fim with resolutions < 2" by Keto et al. (1992) and Campbell et al. (2000). 
Their observations failed to detect component B, suggesting that this source is more deeply 
embedded. They also showed that component C does not have a cometary morphology in 
the mid-IR; the general discrepancy between mid-IR and radio morphologies is discussed 
by Hoare et al. (2007). Campbell et al. (2000) detected two additional sources, which they 
designated E and F. Component F is located about 17" south of component C in the 20.6 //m 
image; component E is located 2" south of and is blended with component C (Campbell et al. 
2000). Mid-IR sources C and E coincide with the arc of the cometary H II region (see Fig. 
9) and were resolved into four individual sources in the sub-arcsecond resolution 10.5 /xm 
and 18.1 iim images of De Buizer et al. (2003). 

The IRAC images of G34.26 illustrate how active and dynamic this MSFR is. Fig. 9 

shows the 4.5 /im and 3.6 fim images. The 8.0 and 5.8 /xm emission is similar to that shown in 
Fig. 9 except that the outflow is only detected at 4.5 fim (see below). The IR counterparts 
of radio components A and C cannot be distinguished from one another, owing to very 
bright, diffuse emission in the IRAC bands. Also, IR components E and C are blended. IR 
component F was detected in all bands. The large ring of emission south-east of components 
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A-C (see left panel of Fig. 9) corresponds to the D radio component of the G34.26 complex. 
Comparison of the 2 cm image of Fey et al. (1994) with the IRAC images shows that this 
IR ring has roughly the same angular extent (~1' radius) and brightness distribution (75% 
complete ring of emission with an opening to the NE) as at radio wavelengths. 

The IRAC images reveal multiple IR sources that were not detected by earher infrared 

observations. Nevertheless, component B remains undetected, even in the IRAC images. 
The point sources listed in the GLIMPSE Catalog and Archive are shown in Fig. 9. The 
high-resolution 2 cm radio contours (Paper I) are overlaid on the 3.6 fim image to indicate the 
positions of the A-C radio components; IR sources E and F arc also marked. Only G34.26 
C has a counterpart in the GLIMPSE Archive (SSTGLMA G034.2572+00.1533) with flux 
densities of 4857 ± 278 mJy and 4481 ± 155 mJy at 4.5 and 5.8 //m, respectively. 

One of the most prominent features in the color composite image of G34.26 (Fig. 9) is an 
outflow detected as an Extended Green Object (EGO, Cyganowski et al. 2009). Molecular 
outflows are particularly strong in the IRAC 4.5 //m band because of shocked H2 and/or CO 
line emission (Cyganowski et al. 2008). The source of the outflow is unknown. Moreover, 
other outflows appear in SiO (Hatchell et al. 2001) and ^^CO (Matthews et al. 1987). The 
latter two outflows are not well-ahgned with the EGO, suggesting that they are different 
outflows arising from different driving sources. 

4.4. G45.07+0.13 

G45.07 consists of a bright NE component (S2cm = 620 mJy) and a much fainter SW 
component (S2cm — 40 mJy), separated by ~6" (see Fig. 6). The SW component was 
flrst detected by Garay et al. (1986) at 2 cm and designated G45.07 A. Both NE and SW 
components are unresolved in our radio continuum maps. High resolution (HPBW ~0('l) 
observations by Turner & Matthews (1984) of the NE component revealed a clumpy ring 
with a gap to the west which they interpreted as a shell with an inner cavity produced by a 
stellar wind from a hot central star. However, Garay et al. (1986) suggested that the ring-like 
structure may be the ionized inner wall of an accretion disk. They further suggest that the 
ring is expanding at ~ 10 km s~^ (owing to a hot stellar wind) and is inclined by ~75° to 
the line of sight. The H76q; RRL integrated over the whole NE component is broad with a 
FWHM of 48.1 ± 0.9 km s"^ (Garay et al. 1986). Garay et al. (1985) also detected broad 
H66q; emission with a FWHM of 42.3 ± 2.3 km s'^. Kcto et al. (2008) reported a detection 
of the H30q; line toward the NE component with FWHM of 33.2 ± 4.2 km s^^ 

We detected H92ci; emission toward both the NE and SW components of G45.07; how- 
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ever, we were not able to determine reliable line parameters and physical properties of the 
sources based on these data. The H92a line observed toward the SW component is very 
narrow relative to the channel width (16.9 km s~^) preventing a reliable Gaussian analysis. 
We estimate that the peak line-to-continuum ratio of this line is > 0.07. For the NE compo- 
nent, the H92a line has a very poor signal-to-noise ratio, likely the result of the high optical 
depth which attenuates the line. The peak line-to-continuum ratio for this line is ~0.03 and 
the deconvolved FWHM of 40 ± 9 km s"^ . 

Hunter et al. (1997) reported a molecular outflow centered on G45.07 NE based on 
observations at millimeter and submillimeter wavelengths and in the CO(6-5) and CS(2-1) 
molecular lines. Infall was also detected by redshifted CS(2-1) absorption, confirming the 
young age of this source. The outflow axis is aligned spatially and kinematically with the 
H2O maser spots in G45.07 NE (Hofner & Churchwell 1996), suggesting that the masers 
form in warm and dense shocked gas at the inner edges of the outflow lobes (Hunter et al. 
1997). 

Several masers have been detected in the NE component of G45.07, while no maser 
emission is associated with the SW component (see Fig. 6b). Hofner & Churchwell (1996) 
detected four clusters of water masers in the area, three of which are closely associated with 
the cm continuum emission from the NE component. The fourth group is located about 2" 
North of the NE component and is not coincident (in high resolution images) with any radio 
continuum emission. The H2O masers are intermixed with 1665 and 1667 MHz OH masers 
at these two locations (Garay et al. 1985; Argon et al. 2000). A 44 GHz GH3OH maser was 
detected ~6" north of the NE component by Kurtz et al. (2004). Single-dish observations 
by Menten (1991) and Caswell et al. (1995a) revealed methanol 6.67 GHz masers. 

4.4- 1- Mid- Infrared Emission 

High-resolution observations (from I'.'l to 1''7; De Buizer et al. 2003, 2005; Krae- 
mer et al. 2003) at mid-IR wavelengths from 10.5 fira to 20.6 /im, revealed three mid-IR 
components associated with the G45.07 MSFR, designated KJK 1, 2, and 3 by Kraemer et al. 
(2003). The brightest component, KJK 1, is coincident with G45.07 NE and contributes 
roughly half of the total (KJK 1-3) flux at 12.5 //m and 20.6 iim (Kraemer et al. 2003). The 
second brightest component, KJK 2, hes 2" north of KJK 1, coincident with the water maser 
group 2" north of the NE component (Hofner & Churchwell 1996, see discussion above); this 
source is a good candidate for a hot molecular core (HMG). Based on the mid-IR luminosities, 
De Buizer et al. (2003) estimated that the putative HMG in G45.07 would contain an 0-type 
star if it is heated by a single star. KJK 2 may drive the outflow reported by Hunter et al. 
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(1997). The third mid-IR component, KJK 3, corresponds to the SW component of G45.07. 

GLIMPSE 3.6 and 4.5 //m IRAC images of the G45.07 MSFR are presented in Fig. 10. 
At IRAC bands, G45.07 is a bright source. The KJK 1 and 2 sources are not well-resolved 
in these images but are clearly apparent in Fig. 10. KJK 3 is very faint at shorter IRAC 
bands, but is bright at 8 /im where PAH emission dominates. 

The GLIMPSE Catalog and Archive sources shown in Fig. 10 reveal the positions 
of (proto)stars that may be responsible for exciting PAH emission. In general there is a 
good correspondence between the distribution of dust and gas in G45.07; however the dust 
emission extends beyond the region of ionized gas emission at 3.6 cm because the PAH 
contribution to this band traces the PDR that surrounds the nebula. The correspondence 
in projected emission does not necessarily imply that ionized gas and dust are well mixed. 
Additionally, the VLA filters out extended emission and much of the extended IR emission 
is instrumental (sidelobes of the Spitzer point spread function), making the interpretation of 
the data more difficult. 



5. Discussion 

5.1. Radio Recombination Lines 

We detected RRLs from all four sources in our sample (see Table 2 for details on the 
observed transitions). Figure 11 shows the relation between the widths (FWHMs) of the 
RRLs and frequency for these sources. Regions classified as HC H ll candidates and UC H ll 
regions are plotted separately in Figure 11a and lib, respectively. The dashed vertical lines 
indicate the rest frequencies of the observed transitions. 

The change of observed Gaussian line widths with frequency (or principal quantum 
number) presented in Figure 11 clearly illustrates the effects of pressure broadening based 
on our data. The limited sensitivity that is typical of very high spatial resolution RRL 
observations rarely permits the Voigt profile to be seen (Roelfsema & Goss 1992). Very 
broad bandwidths and good signal-to-noisc ratios in the weak RRL wings are a prerequiste 
to fit Voigt profiles; our data are not adequate for such a fitting procedure. 

Figure 11a shows that the lines from HC H ll regions (except G34.26 A), particularly the 
lower-frequency transitions, are too broad to be the result of a combination of thermal and 
turbulent broadening only. If an electron temperature of 10,000 K and turbulent velocity 
of 5 km s^^ arc assumed, the contribution from thermal and turbulent motions to the line 
widths would be only ~22 km s^^. The change in line width with transition apparent in 
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Figure 11a strongly suggests that pressure broadening is a dominant broadening mechanism. 
Owing to the rapid change of pressure broadening with principal quantum number n (FWHM 
of pressure broadening oc n^'^; Brocklehurst & Seaton 1972), the lower-frequency lines are 
affected much more strongly by pressure broadening than the higher-frequency lines. At high 
frequencies, pressure broadening becomes negligible, thus line widths should be dominated 
by thermal, turbulent, and bulk motions. Even though pressure broadening depends so 
strongly on n, accurate adherence to the ■n7^'^ dependence is not expected because of other 
dependences on electron density and electron temperature, which vary within Hll regions 
(e.g. Roelfsema & Goss 1992). 

Pressure broadening is less strongly dependent on electron density than on the principal 
quantum number (FWHM oc rie^true, where Ue^tme is the local true electron density; Brockle- 
hurst & Seaton 1972). The true electron density traced by line widths is higher than the rms 
electron density measured from the continuum data. The difference depends on the filling 
factor / of the Hll region (/ = nl^rms/''^e,true)- The filling factor is the fraction of the total 
volume that is occupied by dense gas; thus it is a measure of dumpiness in the H ll region. 
The higher electron density results in more pressure broadening, thus this effect may explain 
the observed increase in line width for the relatively low rms electron density UC H ll regions 
G10.96 E and G28.20 S (see Figure lib). 

Sources G28.20 N, G28.20 S, G34.26 B and G34.26 C have been observed in high fre- 
quency RRLs: H30a for G28.20 N and H53a for the remaining three sources. By comparing 
these high frequency lines with lower frequency (and hence pressure-broadened) lines, we 
can determine ne,true- Comparing these densities with the rms values from the continuum 
observations then provides an estimate of the filling factor. Sources showing the steepest 
slopes in the FWHM versus frequency plot of Figure 11 will have rie^true much greater than 
'>T'e,rms, hcncc havc the smallest filling factors. 

The largest change of linewidth with frequency of the RRL is obtained for G28.20 N. 
Sewilo et al. (2008) obtained an ne,true more than an order of magnitude higher than ne,rms 
and a filling factor of 6 x 10""^ for this source. G28.20 S, on the other hand, shows only a 
4 km s~^ change in linewidth between the H76q; and H53a lines. We determine an Ue^tme 
of (1.2 ± 0.7) X 10^ cm~^ compared to the continuum value of rie^rms = 8 x 10^ cm~^, thus 
yielding a filling factor of 5 x 10~^. 

For G34.26 components B and C we calculate true electron densities of (4.2 ± 0.9) x 10^ 
cm~^ and (4 ± 2) x 10^ cm~^, respectively, resulting in fiUing factors of 0.3 and 0.06. 

We note that the true electron densities and hence the filling factors are calculated under 
the assumption that the H53q; linewidth has no contribution from pressure broadening. If 
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this assumption is not true, then the values we report should be considered as limits; the 
true electron densities will be higher and the true filling factors will be smaller. 

The broad linewidths (FWHM > 40 km s~^) detected toward the HC Hll candidates 
(Figure 11a) probably have a contribution from bulk motions: G28.28 N has been related 
to rotation and possible outfiow; G34.26 B to an outfiow; and G45.07 NE to expansion (see 
Section 4.4). All these phenomena may contribute to line broadening. For the UC Hll region 
034.26 C (Figure lib), Garay et al. (1986) proposed that the relatively broad H76q; line 
results from microturbulent motions of ~30 km s~^ (see Section 4.3). 

5.2. Continuum Emission 

The continuum observations of MSFRs G10.96 and G28.20 at two wavelengths (Ta- 
ble 3) allow a determination of the spectral index for the individual sources in these regions, 
provided that the reliable flux density measurements can be made (GIG. 96 W, G28.20 N). 
For consistency, the two continuum images for each region (1.3 cm and 2 cm for GIG. 96; 7 
mm and 2 cm for G28.2G) were convolved to a common resolution before measuring the flux 
densities. For GIG. 96 W, the 1.3 cm and 2 cm flux densities of 296±15 mJy and 288±15 
mJy, respectively, give a spectral index of -|-(0.07±0.17), indicating that this source is in the 
optically thin regime. The 7 mm and 2 cm flux densities of G28.20 N are 723±37 mJy and 
490±25 mJy, respectively, indicating a spectral index of -|-(0.4±0.1). This value suggests 
that the Hll region has a turnover wavelength in the 2 cm - 7 mm range. The fiux density 
distribution for G28.20 N, based on high-resolution interferometric observations covering the 
wavelength range from 6 cm to 1.3 mm, is shown in Fig. 12. The plot shows the spectrum 
becoming optically thin in the 20 - 50 GHz range. 

We found that G28.2G N is very similar to the well-studied HC Hll region G24.78-F0.08 
Al (Beltran et al. 2007) in terms of physical properties, kinematics (rotation, infall, outflow), 
morphology (shell- like), and association with masers. The flux density distributions of these 
sources are similar in terms of the shape and the turn-over frequency (Beltran et al. 2007; 
Fig. 12). Neither of these HG Hll regions has a rising spectrum to millimeter wavelengths (as 
do G9.62+0.19 E and G75.78+0.34 (H2O), see Franco et al. 2000). The shell-like morphology 
appears to be a common feature of the HG Hll regions. Besides G28.20 N and G24. 78-1-0. 08 
Al, the A and B components in G34.26, and also G45.07 NE, show a shell morphology in 
high resolution images reported by Avalos et al. (2009) and Turner & Matthews (1984), 
respectively. 
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5.3. Distribution of Masers 

We investigated the distribution of maser emission in MSFRs from our sample. We 
found that maser emission is only detected toward sources with broad RRLs (typically the 
most compact and dense sources in the UC/HC Hll region complexes): G10.96 W, G28.20 
N, G34.26 B and C, and G45.07 NE. 

In Figure 4 we show the maser positions for G28.20 N overlaid on the continuum and 
thermal molecular line emission. The OH masers are distributed linearly over a region ~1''4 
(or ~0.04 pc; Figure 4a) with a velocity gradient of ~150 km s"-*^ pc~^. The velocities 
increase from ~94 km at the north-western end of the linear feature to ~100 km s"^ at 
the south-eastern end. The small cluster of maser spots north of the south-eastern end has 
a velocity of ~95 km s~^. The 1.3 cm image (Figure 4c), with 20x higher angular resolution 
than the 7 mm image (Figure 4a), shows the OH masers situated along the south-western 
boundary of the continuum source, at the edge of the CH3CN emission (Figure 4d). The 
location of the OH masers suggests that they arise within a molecular shell surrounding the 
ionized gas. 

The OH maser velocities in G28.20 N roughly agree with the CH3CN velocity of ~95- 
96 km s~^; however, they are redshifted relative to the Hll region velocity of ~93 km s~^ 
estimated by the systemic velocity of the highest frequency RRL observed toward G28.20 
N (H30a at 1.3 mm, Keto et al. 2008; Welch & Marr 1987). A similar redshift is seen for 
CH3OH masers, detected over the velocity range 94 - 105 km s~^. These redshifts for OH 
and CH3OH masers, together with the fact that they appear projected against the optically 
thick continuum, suggest that they arise on the near side of G28.20 N, thus supporting the 
scenario that the masers are located in a collapsing molecular envelope surrounding the H ll 
region. A similar scenario was proposed for W3(0H) (Reid et al. 1980) and several other 
UC Hll regions (including G34.26 C, and G45.07 NE; Garay et al. 1985) based on the maser 
and radio continuum emission in those regions. This interpretation can also apply to G10.96 
W from our sample, where the two CH3OH maser spots are near the continuum emission 
peak, and have velocities (~25 km s~^; Walsh et al. 1998) redshifted by ~14 km s"^ relative 
to the H66a; velocity (see Table 7). 

The OH maser velocities in G28.20 N appear to be inconsistent with the velocity gradient 
reported for the H53a line by Sewilo et al. (2008). In particular, the H53« velocity gradient, 
when extrapolated to the position of the arc of maser emission, gives velocities of 115 - 120 
km s~^, significantly higher than the OH maser velocities of 93 - 101 km s"-*^. The differing 
direction and magnitude of the H53q; and the OH velocity gradients suggests that dynamics 
of the ionized and molecular gas trace different phenomena. 
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The orientation of the H2O masers (SE-NW) in G28.20 N is consistent with the outflow- 
suggested by Sewilo et al. (2008) and Qin et al. (2008). However, the velocity structure of 
the masers is ambiguous and does not provide strong evidence for the outflow. 

In G45.07 NE, the H2O maser spots are ahgned spatiaUy and kinematicaUy with the 
outflow axis detected in CO(6-5) and CS(2-1) molecular hues (Hofner & Churchwell 1996; 
see Section 4.4), suggesting that the masers form in warm and dense shocked gas at the inner 
edges of the outflow lobes (Hunter et al. 1997). 

5.4. Spitzer/IRAC Mid-IR Emission 

We flnd that all the radio components have point source counterparts in the GLIMPSE 
Archive except G34.26 A and B, and G45.07 SW. A faint mid-IR source was detected at the 
position of G45.07 SW; however, it was not included in the GLIMPSE catalog due to the close 
proximity of a much brighter source (G45.07 NE) that prevented an accurate photometric 
measurement. G34.26 A does not have a clear counterpart in the IRAC image; this source 
was detected in the mid-IR by earlier studies, unlike component B which is undetected in 
the mid-IR. 

In Table 10, we list GLIMPSE counterparts to the radio sources (the closest neighbors). 
The distance between the mid-IR and radio positions ranges from 0''2 for G45.07 NE to 2''4 
for G28.20 S (see Table 10). The proximity of the mid-IR sources to the radio peaks of HC 
Hll regions suggests that these sources probably are the ionizing (proto)stars. 

The relatively low density of stars nearby to G10.96 E and G28.20 relative to more 
distant parts of these fields, suggests extinction of background stars in these regions. The 
majority of IR sources shown in Figures 7-10 are probably background stars; however, we 
expect to flnd a number of young stellar objects (YSOs) in these MSFRs. The YSOs can be 
identifled by their mid-IR colors and by fltting their spectral energy distributions with YSO 
models. This would be particularly interesting for G28.20, which is located in an IRDC (see 
Figure 8). 

The Spitzer images of the sources from our sample suggest that generally the large-scale 
spatial distribution of gas and dust in these regions are similar; G34.26 C being the principal 
exception in our sample. Higher resolution infrared observations would be worthwhile, how- 
ever, to separate diff'usc emission from the mid-IR emission from outflows and disks (e.g., 
NGC7538 IRS 1; De Buizer & Minier 2005). 
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6. Summciry 

We present high spatial resolution (0''9 to 2!'3), multi-transition (H53« to H92a) Very 
Large Array observations of four massive star formation regions hosting HC Hll region 
candidates identified in earlier studies. G10.96, G28.20, and G34.26 contain HC Hll region 
candidates with broad H92q; (G10.96 W, G28.20 N) or H76q; (G34.26 B) hnes (Sewilo et al. 
2004a). An additional HC H ll region candidate (G45.07 NE), known to have broad H66q; and 
H76q; lines, small size, high electron density and emission measure, was also included in this 
survey. We observed up to two additional transitions per source to investigate the change 
of line width with frequency. We found a strong correlation of line width with principal 
quantum number, i.e., higher frequency lines of sources with broad H92a or H76a lines 
have narrower linewidths (<40 km s^^) indicating that pressure broadening is the main line 
broadening mechanism at lower frequencies. This result is not surprising because of the 
high electron density in HC H ll regions. However, large scale dynamics is also an important 
contributor to line broadening (Sewilo et al. 2008). We determined physical properties of the 
H II regions based on both the continuum and line data. Multi-frequency RRLs reported in 
this paper are well-separated in principal quantum number and will provide good constraints 
for non-LTE models of these regions. 

We also present a comparison of the radio continuum and mid-IR emission in HC H ll 
regions from our sample. Based on the Spitzer GLIMPSE survey, covering a wavelength range 
from 3.6 to 8.0 /xm, we find a generally good correspondence between the radio continuum 
and dust distributions for the sources in our sample. The mid-IR images, however, show 
more complex features and reveal multiple sources that do not correspond so well to the 
radio continuum. 

We report the distribution of H2O, OH, and CH3OH masers with respect to the radio 
continuum. We also present high-resolution VLA archival data for the water masers in 
G28.20 N for the first time. We compile information from the literature on maser emission 
in G10.96, G28.20, and G45.07. A detailed discussion of the maser distribution in G34.26 can 
be found elsewhere (e.g., Hofner & Churchwell 1996, Gasiprong et al. 2002). As expected, 
all of the massive star formation regions studied in this paper are associated with maser 
emission; however, the masers are not coincident with all the radio components in these 
regions. In G10.96, G28.20, and G45.07 the masers are associated only with the HC H ll 
regions (G10.96 W, G28.20 N, G45.07 NE). In G34.26 the masers are coincident with both 
the HC Hll region G34.26 B and also with the cometary UC Hll region G34.26 C. All of 
the sources associated with masers have broad RRLs. 

The question of the origin of BRRLs is not fully answered, although our observations 
show that pressure broadening is an important mechanism. To distinguish the effects of 
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pressure broadening from large scale motions of the gas and at the same time resolve the 
structures of HC Hll regions will require high-resolution observations of millimeter RRLs 
such as H42q; and H30q; . With its high spatial resolution and sensitivity, ALMA will be the 
instrument to address this issue. 
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Table 1. Instrumental Parameters for the VLA Observations 



Parameter 2003 Mar 20 2004 Mar 7, 8"* 



Program 


AS753 


AS797 


Transitions observed 


H53a 


a) H66a 






b) H76a 






c) H92a 


Observed sources 


G28.20-0.04 


GIO. 96+0.01 (ab) 




G34.26+0.15 


G28.20-0.04 (b) 






G45.07+0.13 (c) 


Total observing time 


10 hr 


a) 3 hr 






b) 9 hr 






c) 1 hr 


Configuration'^ 


D 


C 


Synthesized beam/FWHM ("), PA(°): 






GIO. 96+0.01 




1.4 X 0.8, -8 (a) 






2.0 X 1.2, -2 (b) 


G28.20-0.04 


1.9 X 1.4, -10 


1.5 X 1.4, -10 (b) 


G34.26+0.15 


1.6 X 1.5, -4 




G45.07+0.13 




2.5 X 2.4, -31 (c) 


Rest frequency of the H line 


42,951.97 MHz 


a) 22,364.17 MHz 






b) 14,689.99 MHz 






c) 8,309.38 MHz 


Observing mode 


lA 


2AD 


Bandwidth 


25 MHz 


12.5 MHz 


Number of channels 


31 


31 


Channel separation 


781.25 kHz 


390.625 kHz 


Velocity resolution'^ 


6.5 km s~^ 


a) 6.2 km s""'^ 






b) 9.6 km s'^ 






c) 16.9 km s-i 


Flux density calibrator (Jy): 






3C 286 for GlO.96+0.01 




2.53 (a), 3.50 (b) 


G28.20-0.04 


1.47 


3.50 (b) 


G34.26+0.15 


1.47 




G45.07+0.13 




5.27 (c) 


Phase calibrators (Jy): 






1820-254 for GlO.96+0.01 




0.610 ± 0.004 (a) 






0.692 ± 0.003 (b) 


1851+005 for G28.20-0.04 


0.75 ± 0.02 


0.917 ± 0.003 (b) 


G34.26+0.15 


0.72 ± 0.02 





Table 1 — Continued 



Parameter 


2003 Mar 20 2004 Mar 7, 8^ 


1925+211 for G45.07+0.13 


2.57 ± 0.01 (c) 



^a, b, and c indicate values for the H66a, H76a, and H92a lines, respec- 
tively. 

'^Primary beams are: 1' for AS753; 2', 3', and 5(4 for AS797 (a), (b), and 
(c), respectively. 

'^The observations were done with uniform spectral weighting, thus the 
velocity resolution is 1.2 times the channel separation. 



Table 2. Transitions Observed toward Individual Sources and Adopted Distances 



Source 



H53q; H66a H76a H92q; D (kpc) Ref. 



GlO.96+0.01 
G28.20-0.04 . 
G34.26+0.15 
G45.07+0.13 



+ 



+ 
+ 



+ 

+ 



+ 



14.0 



+1.3 
-0.9 



c, 7 +0.5 

3.7 
6.0 



1 
2 
3 
3 



Note. — The plus signs (+) indicate the observed transitions. 
^Transitions reported in Sewilo et al. (2004) 



References. — (1) Sewilo et al. 2004a; (2) Fish et al. 2003; (3) Araya et al. 
2002 



Table 3. Continuum Parameters^ 



Time Source Integrated 

Position on Diameter Source Flux RMS 



Source A Comp. a(J2000) <5(. 72000) Source 0maj x Smin PA^ Density in Image 

(h m s) (° ' ") (li) (" X ") (°) (mjy) (mjy beam"!) 



GlO.96+0.01 . . . 


1.3 cm 


E° 


18 


09 


42.85 


-19 


26 


29.6 


0.5 












0.4 






W 


18 


09 


39.35 


-19 


26 


28.1 




1.0 


X 


0.9 


177 


291 






2 cm 


E/El'' 


18 


09 


42.98 


-19 


26 


28.2 


1.4 












0.4 






E/E2'' 


18 


09 


42.85 


-19 


26 


30.0 




















W 


18 


09 


39.35 


-19 


26 


28.1 




1.1 


X 


0.9 


53 


289 




G28.20-0.04 .. 


0.7 cm 


N 


18 


42 


58.11 


-04 


13 


57.5 


0.6 


0.7 


X 


0.5 


159 


710 


1.4 






S 


18 


42 


58.12 


-04 


14 


04.9 




2.9 


X 


2.4 


177 


<108<= 








SI 


18 


42 


58.29 


-04 


14 


04.4 












c 






2 cm 


N 


18 


42 


58.11 


-04 


13 


57.4 


2.6 


0.8 


X 


0.6 


150 


494 


0.4 






S 


18 


42 


58.13 


-04 


14 


04.7 




2.9 


X 


2.7 


28 


132 




G34.26+0.15 ... 


0.7 cm 


A 


18 


53 


18.79 


+01 


14 


56.3 


1.0 


0.3 


X 


0.2 


54 


>86<i 


2.8 








18 


53 


18.66 


+01 


15 


00.5 




















Of 


18 


53 


18.56 


+01 


14 


58.2 












6000 s 




G45.07+0.13 ... 


3.6 cm 


NE 


19 


13 


22.08 


+10 


50 


53.2 


0.6 


0.6 


X 


0.6 


87 


350 


0.8 






SW 


19 


13 


21.83 


+10 


50 


48.2 




1.4 


X 


1.1 


82 


60 





''Source positions, deconvolved angular sizes, position angles, and integrated flux densities were determined by elliptical Gaussian fltting 
using the AIPS routine JMFIT. Wc estimate the positional uncertainty to be ~0"1. Uncertainties in source sizes depend on both the brightness 
and the intrinsic source extent; we estimate them to be accurate to within ^10%. The flux density uncertainties are 10% - this value takes 
into account the error of the flt (JMFIT) and systematic effects (e.g. calibration errors, mismatches between the Gaussian model and the 

actual source shape, etc.). 

''The E component of the MSFR G10.96 is overresolved at both 1.3 and 2 cm. Most of the flux is flltered out, thus we do not give either 
the total flux density or source size at these wavelengths. The 1.3 cm coordinates correspond to the maximum pixel. The peak flux density 
at 1.3 cm is 3.9 mjy beam~^. The 2 cm flux density from the region that encloses two compact components (9 a level; see Fig. 2b) is 38 ± 3 
mjy beam"'^. 

■^The faint SI component emerging at 7 mm is partially blended with component S; the reported flux density is the total integrated flux 
density from both components. The flux density of the SI component constitutes ~15% of the total flux. The peak 7 mm flux density of SI 



is 26.6 mjy beam~i. 

'^Thc A component is partially blended with component C; the integrated flux density should be considered the lower limit. 

''The B component is not fully resolved from the cometciry component C, thus we cannot give either its angular diameter or flux density. 
The position given in the table was determined from the 3.6 cm continuum map in Sewiio et al. (2004). 

^The coordinates correspond to the maximum pixel in the cometary component C. The peak flux density is 2.78 Jy beam~^. 

^Thc integrated flux density was calculated using the AIPS task IMEAN from the area of ~22" by ~12" in R.A. and deck, respectively, 
that encloses all three components of the G34.26 H ii region complex. The total flux density of the A component is ~92 mJy, the remaining 
~5908 Jy arise from components B and C. Only a small fraction of this flux density (< 3%) comes from the B component (see text). 
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Table 4. Positions and Observed Parameters of Water Masers in G28.20-0.04 N 



water 


a(J2000) 


(5(J2000) 




C b 

^peak 


maser 


(h m s) 


(° ' ") 


(km s~^) 


(Jy) 


1 


18 42 58.07 


-04 13 56.9 


97.7 


1.2 


2 


18 42 58.12 


-04 13 57.8 


95.6 


10.7 


3 


18 42 58.11 


-04 13 57.6 


93.1 


1.6 


4 


18 42 58.06 


-04 13 57.0 


91.0 


22.2 



^The velocity of the peak channel. 

''The uncertainties of peak flux densities are domi- 
nated by systematic errors. We adopt a standard 10% 
uncertainty for the VLA at K-band. 



Table 5. Source Averaged Physical Parameters based on Continuum Observations^ 



diameter 



Source 




A 


of snhere 


' c 


n<= ™o / 10^ 


EM/10^ 


u 


Mhtt 


loe N' 


Snectral 






(cm) 


fmDcl 




fcm~^) 


(pc cm~"^) 


(pc cm~"^) 






Tvpe 


G10.96+0.0r 


W 


1.3 


121 


0.02 


2.9 


5.3 


57.6 


0.6 


48.8 


07.5V 






2 


124 


0.06 


2.7 


5.3 


56.4 


0.6 


48.8 




G28.20-0.04 


N 


0.7 


28 


0.04 


17 


36 


43.9 


0.05 


48.4 


08V 






2 


35 


0.2 


9.8 


19 


37.1 


0.05 


48.2 






S 


0.7 


131 


0.0003 


0.7 


0.3 


23.4 


0.2 


47.6 


09.5V 






2 


120 


0.01 


0.8 


0.4 


23.9 


0.2 


47.6 




G34.26+0.15'^ 


A 


0.7 


8 


0.03 


27 


24 


16.3 


0.002 


47.1 


B0.5V 


G45.07+0.13 


NE 


3.6 


32 


0.7 


9.4 


19 


33.4 


0.04 


48.1 


08.5V 




SW 


3.6 


65 


0.03 


1.4 


0.8 


18.7 


0.05 


47.3 


BOV 



Note. — The physical parameters are: the linear diameter of a spherical Hll region (assumed model), the continuum 
optical depth (tc), the rms electron density (ne^rms), the emission measure (EM), the excitation parameter (U), the 
mass of ionized gas (Mhh), and the Lyman continuum photon flux required to maintain ionization of the nebula (N^). 

'^We estimate an uncertainty of the derived physical parameters to be up to 20%. These parameters depend on the 
different combinations of the integrated flux density, observed source size, and assumed distance. All of these quantities 
have uncertainties of ~10%. 



'^Lower limits to the spectral types of the exciting stars are derived from using the model stellar atmosphere results 
of Smith et al. (2002). See text for details. 

^We do not report the physical parameters for component E of GlO.96+0.01, because most of the emission from this 
source if filtered out. See footnotes to Table 3. 

'^We do not report the physical parameters for components B and C of G34.26+0.15. They are blended at our spatial 
resolution, hence we cannot determine their angular sizes or flux densities. 



Table 6. Peak Values of the Physical Parameters^ 



Source 




u 




As 


Tb 




EM/10^ 


■^e,rms/lO 






(GHz) 


(mJy/beam) 


(pc) 


(K) 




(pc cm~^) 


(cm~^) 


VjriU.yb+U.Ui 


W 


ZZA 


1 /1 7 o _i_ n /I 
14/. 2 ± U.4 


0.0 / 


OO 

oZo 


n no 
O.Oo 


7 1 
I.L 


o o 
o.o 






1/17 


loo.O ± U.4 


U.U/ 


4oU 


n n 
U.UO 


A O 

4.Z 


O A 

ZA 




Hi 


ZZA 


d.y ± u.4 




n 

y 


\n nni 
>U.UU1 


n o 
U.z 






T71 b 


1 A '7 

14. / 


D.O ± U.4 




lb 


n nno 
U.UUz 


n 1 
U.l 








1 /I "7 

14. / 


D.y ± U.4 




1 / 


n nno 
U.UUz 


n 1 
U.l 




r^OQ on n n/i 
kjiZo.Zu — U.U4 


1\T 


AO n 
4z.y 


DOl ± 1 


U.Ulo 


IDO 


n ni 7 
U.Ul / 


1/1/1 
14.4 


n Q 

y.o 






14.7 


401.9 ± 0.4 


0.02 


1079 


0.1 


10.1 


7.3 




s 


42.9 


27 ± 1 


0.07 


7 


>0.001 


0.6 


0.9 






14.7 


30.7 ± 0.4 


0.08 


82 


0.01 


0.7 


1.0 


G34.26+0.15 


A 


42.9 


86 ± 3 


0.004 


24 


0.002 


2.1 


7.1 




B 


42.9 


120" ± 3 




33 


0.003 


2.9 






C 


42.9 


2777 ± 3 




773 


0.1 


69.7 




G45.07+0.13 


NE 


8.3 


328.3 ± 0.7 


0.02 


988 


0.1 


2.7 


4.0 




SW 


8.3 


49.4 ± 0.7 


0.04 


149 


0.02 


0.4 


1.1 



^We estimate an uncertainty of the derived physical parameters to be up to 20%. These param- 
eters depend on combinations of the integrated flux density, observed source size, and assumed 
distance. All of these quantities have uncertainties of order 10%. 

'^Two maxima/subcomponents can be distinguished in the 2 cm image of component E of 
GlO.96+0.01 (see Fig. 2). El and E2 correspond to the NE and SW maxima, respectively. 



'^The peak flux density of G34.26+0.15 component B is an upper limit because of blending with 
component C. 



Table 7. Gaussian Fits to the H53q;, H66q;, H76q;, and H92q; Radio Recombination Lines 



Source RRL Vlsr FWHM^ Sl Sc |^ / |^ 









(km s 


-1 \ 

) 


1 


(km s" 


-1 


) 


(mJy) 


(mJy) 




(peak) 




(km s ) 


GlO.96+0.01 . 


E 




11.4 


± 


0. 


,7 


18.7 


± 


1. 


,8 


12.6 ± 1.0 


46.3 ± 2.3 


0. 


,27 ± 


0. 


,03 


5.4 ± 0.7 








15.8 


± 




n 


25.3 


± 


2, 




9 7 + 08 


53 8 + 2 7 

'.J • ) . _1_ ^ . i 


n 


,18 ± 


0. 


02 


4 9 + 07 




W 

V V 


J- J-VJVJCt. 


13.6 


± 


n 


A 


28.5 


± 


1 


n 


359 + 1 1 


147 _)_ 7 


n 


.24 ± 


0, 


01 


7 4 + 05 

1 . t: _1_ \J . tj 






H76a 


13.3 


± 


0. 


,5 


29.6 


± 


1, 


.3 


19.9 ± 0.7 


131 ± 7 


0, 


.15 ± 


0, 


.01 


4.8 ± 0.4 


G28.20-0.04 . 


N 


H53q; 


88.5 


± 


0. 


,5 


39.0 


± 


1 


.4 


144 ± 4 


372 ± 19 


0, 


.39 ± 


0, 


.02 


16.1 ± 1.1 






R76a 


80.1 


± 


0. 


,8 


55.7 


± 


2, 


.2 


17.0 ± 0.5 


235 ± 12 


0, 


.07 ± 


0, 


.01 


4.3 ± 0.3 




s 


H53q; 


99.6 


± 


0. 


,4 


19.9 


± 


0, 


.9 


63.6 ± 2.4 


65.3 ± 3.3 


0, 


.97 ± 


0, 


.06 


20.7 ± 1.6 






H76q; 


99.3 


± 


0. 


,5 


24.8 


± 


1 


.3 


7.3 ± 0.3 


34.9 ± 1.7 


0, 


.21 ± 


0, 


.01 


5.6 ± 0.5 


G34.26+0.15 . 


A 


H53a 


50.2 


± 


0. 


,4 


21.9 


± 


1, 


.0 


46.3 ± 1.7 


62.9 ± 3.1 


0, 


.74 ± 


0, 


.05 


17.2 ± 1.3 




B 


H53a'= 


64.9 


± 


1. 


,0 


34.8 


± 


2, 


.6 


20.5 ± 1.2 


58.6 ± 2.9 


0, 


.35 ± 


0, 


.03 


12.9 ± 1.4 




C 


H53a 


47.9 


± 


1. 


,2 


49.4 


± 


3, 


.7 


589 ± 32 


2587 ± 129 


0, 


.23 ± 


0, 


.02 


12.0 ± 1.2 



'^FuU width at half-maximum. 

^Sc is a continuum flux density of the same area of the source over which a RRL was integrated. It was 
determined using AIPS task IMEAN. 

'^The spectrum was Manning smoothed. 
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Table 8. Corrected H53q;, H66q;, H76q;, and H92a Radio Recombination Line Parameters 

and the LTE Electron Temperatures 



Source 




RRL 


a 

(km s""*") 


o b 
(Sl) 
\Sc ^corr 

( neak 1 


rp* c 

(K) 


GlO.96+0.01 . 


E 


H66a 


15.5 ± 2.1 




11400 ± 1600 






H76q; 


19.6 ± 3.2 




8800 ± 1400 




W 


H66q; 


27.8 ± 1.1 


0.25 ± 0.01 


7500 ± 500 






R76a 


28.0 ± 1.4 


0.16 ± 0.01 


7400 ± 500 


G28.20-0.04 . 


N 


H53q; 


38.5^ ± 1.5 


0.39 ± 0.02 


7000 ± 400 






H76a 


54.9"^ ± 2.3 


0.08 ± 0.01 


7300 ± 500 




S 


H53a 


18.8 ± 0.9 




5900 d= 400 






H76a 


22.9 ± 1.4 




6800 ± 500 


G34.26+0.15 . 


A 


H53a 


20.9 ± 1.0 


0.75 ± 0.05 


6800 ± 500 




B 


H53q; 


33.1 ± 2.8 




8800 ± 900 




C 


H53q; 


49.0 ± 3.7 


0.24^ ± 0.02 


8700 ± 800 



^FuU width at half-maximum (FWHM) of the spectral lines corrected for 
broadening introduced by spectrometer and smoothing. The convolving func- 
tion in the spectral frequency domain was assumed to be a Gaussian function 
with FWHM of 1.2 times a channel width (or twice a channel width if a spec- 
trum is Hanning smoothed) and deconvolved from the observed line widths. 
Velocity resolutions for all frequencies are given in Table 1. 

''We adopt the same uncertainties as for the fitted line-to-continuum ratios 
given in Table 7. 

^For optically thin sources, the values of (^)peak from Table 7 were used to 
calculate T*. 

^For G28.20 N, the H53q; and H76q; hne widths are shghtly smaller that 
those reported in Sewilo et al. (2008) based on the same data: 39.7 ± 1.3 
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km and 57.6 ± 2.2 km for H53q; and H76q; lines, respectively. The 
difference in linewidths is due to different integration areas used in both papers. 
While Sewilo et al. (2008) integrated the lines over the entire source, the lines 
reported in this paper were integrated over the area enclosed by the 50% 
contour. 

''The line-to-continuum ratios corrected for opacity for this source was de- 
rived from the peak optical depths listed in Table 6. 
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Table 9. Flux densities of G28.20-0.04 N 



1/ 


Iiistr. 


kj y xj.ijxj.« 1 i\ ,f III 1 1 


S 


Ref. 


(GHz) 










5.0 


VLA 


2.0x1.6 


150±15 


1 


6.7 


ATCA 


~1.9 


326^ 


2 


8.3 


VLA 


0.9x0.8 


297±45 


3 


8.3 


VLA 


4.9x2.7 


300±30 


4 


14.7 


VLA 


1.5x1.4 


494±50 


9 


14.7 


VLA 


0.6x0.5 


543±81 


3 


22.4 


VLA 


0.3x0.2 




5 


22.4 


VLA 


0.1x0.1 


630±63 


9 


43.0 


VLA 


1.9x1.4 


710±70 


9 


43.0 


VLA 


0.17x0.13 


645±65 


6 


231.9 


SMA 


1 


720'^'^ 


7 


231.9 


SMA 


0.6x0.2 


890±30^ 


8 



'^Uncertainties are not provided in the referenced 
papers. 

^The 1.3 cm flux density given by Sollins et al. 
(2005) seems to be anomously high relative to other 
measurements (see Fig. 12). 

■^Both 231.9 GHz flux densities are from the same 
SMA observations. Any difference in flux density 
represents a difference in the data reduction. 

References. — (1) Purcell ct al. 2008; (2) 
Walsh et al. 1998 ; (3) Kurtz et al. 1994; (4) 
Sewilo et al. 2004a; (5) Sollins ct al. 2005; (6) 
Sewilo et al. 2008; (7) Keto et al. 2008; (8) Qin et al. 
2008; (9) this paper 



Table 10. GLIMPSE Counterparts to the Radio Continuum Sources 



Source 


GLIMPSE ID^ 








Flux densities^ (mJy) 




y ) 


(mpc) 




3.6//m 


4.5//m 


5.8//m 


8 C\iim 


G10.96 W 


SSTGLMA GOlO.9584+00.0219 


2 


136 




9.7 ± 0.9 


59 ± 8 


174 ± 12 


395 ± 31 


G10.96 E2 


SSTGLMA GOlO.9646+00.0098 


0.3 


20 




8 ± 3 


13 ± 4 


59 ± 8 




G28.20 N 


SSTGLMA G028.2003-00.0493 


2.2 


61 


9.0 ± 0.4 


358 ± 19 


2123 ± 99 


3892 ± 99 




G28.20 S 


SSTGLMAG028. 1984-00.0501 


2.4 


66 




2.1 ± 0.3 


9 ± 1 






G34.26 C 


SSTGLMA G034.2572+00.1533 


0.7 


13 






4857 ± 278 


4481 ± 155 




G45.07 NE 


SSTGLMA G045.0712+00.1321 


0.2 


6 


40 ± 1 


2240 ± 184 









^The Galactic coordinates of the Spitzer sources are part of the GLIMPSE IDs. ' 

CO 

•^The distance between the radio continuum source and the nearest mid-IR GLIMPSE point source; kinematic distanbes to 
MSFRs hsted in Table 2 were used to calculate distances in mpc. 

^None of the sources has J and H fluxes available in the 2MASS Point Source Catalog. 
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0.15 
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0.05 
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Fig. 1. — The 1.3 (a) and 2 cm (b) VLA continuum images and the H66q; (c) and H76q; (d) 
Unc profiles integrated over the central portion of the W continuum component of the MSFR 
G10.96. The positions of the two CH3OH maser spots in G10.96 W are superimposed on 
the 2 cm image (b). Beam sizes are the same as in Figure 2. 
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Fig. 2. — The 1.3 (a) and 2 cm (b) VLA continuum images and tlie H66q; (c) and H76a 
(d) line profiles from the E continuum component of the MSFR G10.96. The beam sizes 
shown in the lower left corners of the continuum images are and ~1''6 at 1.3 and 2 cm, 
respectively (see Table 1). The 1.3 cm peak of maximum emission corresponds to one of the 
two peaks seen in the 2 cm image (E2 at south-west, see Table 6). 
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Fig. 3. — (a) The 7 mm VLA continuum image and the integrated H53a hne profiles from 
the N (c) and S (e) components of MSFR G28.20. (b) The 2 cm VLA continuum image and 
the integrated H76q; hne profiles from G28.20 N (d) and S (f). The beam sizes shown in the 
lower left corners in (a) and (b) are ~l'/6 and ~1''4, respectively (see Table 1). 
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Q28.20-0.04 - OH, H2O, and CH3OH maser emission 
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Fig. 4. — The distribution of OH and CH3OH masers in G28.20 N superimposed on: the 7 mm 
VLA image (a and b: b is a magnification of the central portion of the image shown in a; this 
paper), the 1.3 cm VLA image (c; VLA archival data, project AZ168), and the CH3CN SMA image 
(d; Qin et al. 2008). The synthesized beams of the VLA shown in the lower left insets in (a) and 
(c) are ~1{'6 and 0''09, respectively (see Table 1). The synthesized beam of the SMA is ~1''5 (see 
Section 4.2) and is shown in the lower left corner in (d). 
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Fig. 5. — The 7 mm VLA continuum image (a) and the integrated H53q; hne profiles from 
the A (c), B (b), and C (d) components of the MSFR G34.26. The beam size shown in the 
lower left corner in (a) is ~1('5 (see Table 1). 
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Fig. 6. — (a) The 3.6 cm VLA continuum image of the NE and SW components of the 
MSFR G45.07. (b) The magnification of (a) centered on G45.07 NE. Locations of 22 GHz 
H2O, 44 GHz CH3OH, 1665 MHz and 1667 MHz OH masers in the vicinity of G45.07 NE 
are indicated. The VLA beam is shown in both images (HPBW~2''5, see Table 1). 



Fig. 7. — The IRAC 8.0 {left) and 3.6 /iin (right) images of G10.96. The positions of radio 
sources are indicated with red arrows in the 8.0 fim image, while the 3.6 cm radio contours 
are overlaid on the IRAC 3.6 /xm image. The images are centered on the same coordinates 
and have the same size. The sources from the GLIMPSE Catalog and Archive are marked 
in the 3.6 /im image with green and red circles, respectively. The intensity ranges from 0.5 
mJy to 29.2 mJy for the 8 fim image, and from ~1 fiJj to 22.4 mJy for the 3.6 /^m image. 
Thirty arcesc are equal to ~2 pc assuming a distance of 14 kpc (see Table 2). 
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Fig. 8. — The IRAC 8.0 fim (left) and 3.6 /im (right) images of G28.20. The 8.0 /iin image 
shows that G28.20 is located toward an infrared dark cloud (see text). The 3.6 /im image 
is zoomed in on the source; 7 mm radio contours indicate the positions of radio sources 
(see Figure 3a). The point sources from GLIMPSE are marked. The source appearing to be 
blended with G28.20 N slightly to the south in the 8 fim image is an image artifact produced 
by the "electronic bandwidth effect" (see Section 2.2). The intensity ranges from 0.5 mJy 
to 51.2 mJy for the 8 /xm image, and from ~7 /iJy to 25.4 mJy for the 3.6 /zm image. One 
arcmin and 10 arcsec are equal to ~1.7 pc and ~0.3 pc, respectively, assuming a distance of 
5.7 kpc (Table 2). 
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Fig. 9. — Left: The three color composite image of G34.26. Red, green, and blue correspond 
to IRAC 8.0 /im, 4.5 fim, and 3.6 /im, respectively. The 7 mm 3a radio contour (see Figure 
5) is overlaid to indicate the radio positions. A striking feature of the image is a massive 
outflow detected at 4.5 /im (see Section 4.3.2). An incomplete ring of emission dominating 
the central and lower parts of the image corresponds to radio component D. Right: The IRAC 
3.6 /im image of G34.26 with the high-resolution 2 cm radio contours (Sewilo et al. 2004a) 
overlaid. Mid-IR components E and F from Campbell et al. (2000) and the peak of radio D 
component (Fey et al. 1994) are marked with blue circles and labeled. The GLIMPSE point 
sources are marked. The 8.0 fim image suffers from the "electronic bandwidth effect" (see 
Section 2.2). The intensity ranges are: (~1 /iJy, 20.7 mJy), (~1 yuJy, 23.4 mJy), and (0.3 
mJy, 55.9 mJy) for 3.6, 4.5, and 8.0 fim image, respectively. One arcmin and 10 arcsec are 
equal to ~1.1 pc and ~0.2 pc, respectively, assuming a distance of 3.7 kpc (see Table 2). 
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Fig. 10. — The IRAC 4.5 fim {left) and 3.6 fim (right) images of G45.07. The images have 
the same size. The 3.6 cm radio contours are overlaid on both images and the NE/SW 
components are indicated with red arrows in the 4.5 /zm image. Sources from the GLIMPSE 
Catalog and Archive are marked in the 3.6 fim image with green and red circles, respectively. 
The intensity ranges from ~1 /xJy to 28.8 mJy for the 4.5 fim image, and from ~1 /iJy to 
21.6 mJy for the 3.6 fim image. Ten arcsec are equal to 0.3 pc at a distance of 6 kpc (see 
Table 2). 
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Fig. 11. — Plots showing the relation between the half-power widths of the RRLs (FWHMs) 
and frequency for HC H ll region candidates (left) and UC H ll regions (right) from our sample 
(Paper I and this paper). The plots include only the sources observed in more than one RRL. 
The FWHMs with the uncertainties, as well as references can be found in Section 4. The 
values measured in this paper are also listed in Table 8. See the discussion of the plot in 
text (Section 5.1). 
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Fig. 12.— Spectral energy distribution of G28.20 N from 5 GHz to 231.9 GHz. Flux densities 
and beam sizes, along with references, are listed in Table 9. The uncertainties are plotted 
when available. 



